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I. INTRODUCTION 



DLCQ (Discretized Light-Cone Quantization) [|T]-Q| is a suggested computational pro- 
cedure in which one specifies quantization conditions on the characteristic surface x + = 
(x° + x 1 ) = 0, introduces periodicity conditions to induce a discrete basis, truncates the 
basis set by some procedure to produce a finite matrix, then takes the spectrum and eigen- 
vectors of that matrix as an approximation to the physical spectrum and wave functions. 
The difficulty is that, as always in quantum field theory, removing an infinite set of high 
(bare) energy states induces a renormalization of the operators. A completely consistent pro- 
cedure for performing the truncation and renormalization has not yet been demonstrated, 
but the problem has received considerable study. Some calculations have been published Q 
which simply use the periodicity conditions combined with a momentum cutoff. While the 
numerical results are accurate for superrenormalizable theories such as 1 + 1-dimensional 
gauge theories, it is clear that such a procedure is problematical for renormalizable theories 
in 3 + 1 dimensions. A systematic renormalization procedure has been proposed by Wilson, 
Perry and co-workers ||. These authors use a cutoff chosen for computational convenience 
and then try to find the mixing of the operators under renormalization using ideas of the 
Wilson renormalization group. Since the procedure makes no attempt to preserve the sym- 
metries of the theory, one expects a great deal of mixing and many counterterms; some skill 
in guessing the counterterms seems to be required. 

In this paper we will suggest a procedure which lies somewhat closer to traditional ideas 
in field theory than the Wilson plan in that we will make an attempt to preserve more of 
the symmetries of the theory; but the objective is the same. The idea is to add enough 
Pauli-Villars fields |J to the theory to regulate perturbation theory. Having done that we 
hope that, since the theory is basically finite, the periodicity conditions and momentum 
cutoff will be sufficiently benign to allow a consistent renormalization to be performed. In 
the Wilson language, we hope that the heavy fields will add the necessary counterterms 
automatically with no particular cleverness from us. 

In the next section we consider the one- loop fermion self mass in Yukawa theory |7j. This 
problem has been considered previously in the literature ||, but we shall discuss the analysis 
in the context of using the Pauli-Villars program to preserve the discrete chiral symmetry 
of the theory. The computation requires three Pauli-Villars fields for proper regularization, 
including the elimination of all terms — including spurious finite terms — not proportional 
to the bare fermion mass squared. 

We then present and solve a model field theory very similar to the scalar field model 
studied in the 1950's by Greenberg and Schweber |J. This model, which requires renor- 
malization, allows us to illustrate the procedure and to examine some important numerical 
issues, at least within the context of the model. These issues include the number of states 
which must be devoted to the heavy fields and the related question of how heavy their masses 
must be. Section [TV] discusses the numerical solution of this same model and compares these 
results with the analytic solution. 

A final section contains our general conclusions. This is followed by two appendices that 
provide details of our light-cone conventions and of improved methods for accurate DLCQ 
calculations. 
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II. REGULARIZATION OF THE FERMION SELF-ENERGY IN LIGHT-CONE 

QUANTIZATION 



A. Analysis 

We consider Yukawa theory defined by the action 



S 



d x 



(2.1) 



For the problem of interest here, the A0 4 interaction will not be needed. The operator P 
which controls the dynamics is 



P~ 



dx d x T 



(2.2) 



where 



+2^°(-iYdi + M + g<p)i) + + h.c. 



The field is nondynamical and can be eliminated via the constraint relation 

1 



id-ip- 



;7 



-ifdi + M + g<j>il>+. 



(2.3) 



(2.4) 



For the second order shift in the eigenvalue of the operator P of the one-fermion state with 
p± = 0, one easily calculates 



a 

2^2 



f 1 dx r 2 qj + y-xrM' 
Jo 1 - x J q\ + x 2 M 2 + (1 - x)fi 



q 2 ± + (2- xfM 2 



where 



a = 



An 



(2.5) 



(2.6) 



The integral is divergent in the ultraviolet and must be regulated. Let us first consider 
regulating the integral with a momentum cutoff. While several possibilities might be consid- 
ered, including a cutoff on q± alone, the most commonly used cutoffs couple q + and q± in 
some way. To retain boost invariance we will consider the "invariant mass" cutoff in which 
the total invariant mass of the intermediate state is limited Hi~0|l . For the present case this 
rule gives 



ql + V 2 ql + M 2 



+ 



x 



1 



< A 2 



x 



(2.7) 



This cutoff also appears if we simply limit the change in mass of the matrix elements of the 
interaction Hamiltonian. For the integral we then get 



3 



/(/z 2 ,M 2 ,A 2 ) EE 



a r L + dx 



ql + (2- *) 2 M 2 



d q± 



2tt 2 Jl- 1 - x J q l<L q\ + x 2 M 2 + (1 - x)/i 2 



(2.8) 



where 



2A 2 



A + /r — M ± J (A 2 + /x 2 - M 2 ) 2 - 4AV 



L = A 2 x(l - s) - /i 2 (l - x) - M 2 x . 



(2.9) 



The integral can be done in closed form, but the result for arbitrary parameters is not 
terribly illuminating. If we take A 2 ^> /i 2 ^> M 2 we get 



a 



/(/z 2 ,M 2 ,A 2 )^- 
In 



y-/i 2 lnA 2 + /i 2 ln/i 2 -^ 



+M 2 3 In A 2 - 31n/x 



2 9 5/x 



2 



M 4 (^ln(M 2 /^ 2 ) + — - 




(2.10) 



Perhaps the most striking thing about this result is not so much that it is quadratically 
divergent as A — > oo, but that it does not go to zero with M. This is in contrast to 
the Feynman result. In fact the vanishing of the self mass with vanishing bare mass is 
formally protected by the discrete chiral symmetry: ip —>■ ij^ip, with <p — ► —cf). That I is not 
proportional to M 2 is due to the fact that the invariant mass regulator does not preserve 
this symmetry. 

It is well known that the four-dimensional Feynman integral for the fermion self energy 
can be regulated by the addition of one Pauli-Villars boson field. If that is done the integral 
is then "finite" by power counting and vanishes with M. One might then think that if one 
first performs the q~ integral one would get a finite three-dimensional light-cone integral. 
However, the Feynman integral is only conditionally convergent, and therefore any value 
ascribed to it is a prescription. The standard integration procedure - symmetric integration 
in the spatial momenta with the q° integral done last - preserves the discrete chiral symmetry 
and thus leads to the vanishing of the result at M = 0. The terms in ( 2.10Q which do not 
vanish as M — > at large A include terms quadratic in A, logarithmic in A, and independent 
of A. Therefore three Pauli-Villars bosonic fields are necessary to render the light-cone 
integral consistent with discrete chiral symmetry. The entire light-cone integral (|2.10|) is 
then finite and vanishes as M — > 0. Thus we need three Pauli-Villars conditions: 



(2.11) 



where the aij's and /Vs are the coupling constants and masses of the heavy fields. The 
logarithmic divergent term 3M 2 In A 2 in ( [2.10|) returns if the masses of the heavy fields go to 
infinity, but the nonzero value at M = does not return. The fact that three Pauli-Villars 
fields are necessary to regulate the self-energy graph in the light-cone representation is an 
old result [O, and it has received considerable study in [T3]. One might wonder whether 
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there is some feature of the theory from a purely equal-time perspective that would allow one 
to predict the number of heavy fields necessary to regulate the calculation in the light-cone 
representation. [] 

To perform a DLCQ calculation one must limit the range of the momenta (which the 
A cutoff does) as well as make the momenta discrete by introducing periodicity conditions 
for the fields on the surface x + = 0. We may take ifi + = A + ip to be antiperiodic in x~ and 
periodic in x±. (See Appendix |A|.) The scalar fields are taken to be periodic in both x~ 
and Xj_. With the heavy fields and the momentum cutoff in place, the only effect of the 
periodicity conditions on the above perturbation calculation is that the finite integral is then 
evaluated as a discrete sum. The convergence of the discrete sum to the continuum result 
is discussed below. 

The discussion in this section suggests a new general procedure for resolution of the 
UV divergences of light-cone Hamiltonian field theory. The heavy fields will produce the 
counterterms necessary to make a consistent renormalization possible. What we propose is to 
test this procedure nonperturbatively, that is, include enough heavy fields in the Lagrangian 
to regulate perturbation theory, then produce a finite matrix with a momentum cutoff 
and discretization. The periodicity conditions can also lead to constrained zero modes, as 
discussed in [13]. Another important advantage of the Pauli-Villars fields is that the terms 
from the constrained zero modes which would affect the one loop mass shift (the most 
singular terms due to the constrained zero modes) are zero at the level of P~. In the present 
paper we shall not attempt a full Yukawa calculation. We shall illustrate the method for a 
soluble model and provide some discussion of numerical issues. 



B. Discrete Evaluation of the Light-Cone Integral 

1. DLCQ 

From the periodicity conditions in the light-cone box 

-L<x~<L, -L ± <x,y<L±, (2.12) 
one obtains discrete momenta 

i 7T .7T7T. . . 

p + ^-n, p ± -> (— n x , — n y ) , (2.13) 

Li ±jj_ ±j_\_ 

with n even for bosons and odd for fermions. Integrals are then replaced by discrete sums 
obtained as trapezoidal approximations on the grid of momentum values 



|^ + /rfW(p + ,P±)^^(^)'E E f(nn/L,n ± n/L ± ). (2.14) 



2 N_ 



x We are not certain, but we speculate that it may be the same number necessary to regulate 
time-ordered perturbation theory in the equal-time representation; perhaps the results of on 
the connection between the two representations could then be extended to higher order. 
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The limit L — ► oo can be exchanged for a limit in terms of the integer resolution |l] K = ^-P + . 
The longitudinal momentum fraction x = p + / P + becomes n/K. H^c is independent of L. 

Because the longitudinal integers n are always positive, DLCQ automatically limits the 
number of particles to no more than ~K/2. The integers n x and n y range between limits 
associated with some maximum integer N±_ fixed by L±_ and a cutoff that limits transverse 
momentum. 

The momentum-space continuum limit is reached when K and N± become infinite. The 
transverse length scale L± is chosen such that N±n / L± is the largest transverse momentum 
allowed by the cutoff. The integrations for Pauli-Villars subtractions use the transverse 
scale L± determined for the physical boson. This insures use of a common grid that can 
easily represent momentum conservation in interactions. 

We then compute the dimensionless integral 

/(/?, M 2 , A 2 ) = —J(fi 2 , M 2 , A 2 ) (2.15) 

and the subtracted integral 

WM 2 /^ 2 , ^//i 2 , A 2 //, 2 ) ee IV, M 2 , A 2 ) + £ M\ A 2 ) . (2.16) 

i a 

When the loop integrals are evaluated numerically, the proportion of error increases with 
each subtraction. Once all three Pauli-Villars subtractions are done, the error can dominate. 
This is the case for an ordinary DLCQ calculation. The individual integrals are large and 
therefore must be computed accurately for the differences to be accurate. It can be helpful 
to have well separated Pauli-Villars masses, because the coefficients oti are then 0(1) and 
do not amplify the errors. 

The domain of integration, as defined by the invariant-mass cutoff, is not commensurate 
with the DLCQ grid. This causes errors of two types: one is a truncation error where the 
edge of the domain is not properly counted and the other is the loss of rotational symmetry in 
the transverse grid. In one-dimensional cases, only the former type occurs. In that context, 
it is easily handled as part of an extrapolation in K; in the context of the three iterated 
integrals used for three dimensions, the error becomes much less controlled. 

If the trapezoidal rule is applied to each iterated integral, as is done in the standard 
DLCQ approach, the errors will not follow a systematic dependence on the grid spacings for 
a reasonable number of grid points. This lack of systematic dependence on K and N± can 
be seen in Fig. [I]. 

These errors could be overcome with a commensurate grid that uses polar coordinates 
in the transverse direction. This would not be easily extended to situations with more than 
two particles. Also, it turns out that although a commensurate grid controls the errors in a 
systematic way, the errors are still large. Other methods that use the DLCQ grid have been 
found superior. 

Given the rectangular DLCQ grid, one can improve on the simple application of the 
trapezoidal rule used in (|2.14j ). The alternative integration schemes are of the general 
form 

Jd n rf(r)^Y.^,U(r>, J ,...), (2.17) 
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145 




FIG. 1. One- loop fermion self energy. The horizontal line is the exact result. The smoothly 
curved results come from use of transverse circular weighting and longitudinal Simpson weighting. 
The scattered results are from ordinary DLCQ calculations. The DLCQ grid parameters take the 
ranges K = 10, 12, . . . , 24 and N± = 5, 6, . . . , 30. The lines connect points calculated with the 
same K value. 
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TABLE I. Values of the subtracted integral I su b(M 2 //i 2 , n 2 /^ 2 , A 2 /^ 2 ) in the limit of infinite 
cutoff. The Pauli-Villars masses are [if = 10/i 2 , /i| = 50^i 2 and /j, 2 = WO/j, 2 . 





M 2 : 





O.OLu 2 


0.05/i 2 


O.lfj? 


0.2/i 2 


DLCQ, 
Exact, 


improved and extrapolated: 
to order M 4 and M 4 lnM 2 : 


-0.064 
0.0 


0.11 
0.1402 


0.70 
0.6661 


1.37 
1.2721 


2.70 
2.3778 



where, unlike the case of the trapezoidal rule, the weights ifjj,... will not all be equal. Special 
formulas for intervals near the edge can be derived, and one can even consider variations 
on the higher-order Simpson's rule. The transverse integrations can be treated in polar 
coordinates with a basis of circles of irregular radii chosen to pass through the points of the 
square grid. These methods provide results for integrals far better than ordinary DLCQ, as 
shown in Fig. [I]. A discussion of the details is given in Appendix [B|. 

2. Results 

The results for the momentum-space continuum limit of the discrete sums are obtained 
by extrapolations which use values of 20, 22, and 24 for K and 25 through 30 for N±. All 
results are given in units of the boson mass fi. These are fit by least squares to either 
c + cii/K 3 + h/Nl or c + ai/K 3 + a 2 /K A + b x /N\ + b 2 /N\. The latter is used for the /ii 
integral. This means that at most 5 parameters are used in fits to 18 points. 

Extrapolation to the continuum after subtraction is not as accurate as extrapolation of 
each integral separately. The subtraction of the discrete sums induces a greater variation in 
errors that is harder to fit properly. 

The range in N±_ was selected to avoid values where the /12 integral was badly approxi- 
mated. However, here "badly" is to be interpreted relative to the desired final error of 0.02, 
which is slightly more that 0.2% of the answer. 

Values of the subtracted integral for different fermion masses M are plotted as functions 
of 1/A 2 in Fig. |2|. The extrapolation to A = 00 can be done by fits to 1^ + a/ A 2 . They yield 
the values in Table |I[ In obtaining these values, the error in each individual integral has been 
reduced to ±0.02 as measured against the analytic result at M 2 = 0. This implies an error of 
±0.04 in the subtracted result. Extrapolation in A 2 induces additional uncertainty reflected 
in the miss by 0.06 of zero for M 2 = 0. The ratios of the tabulated values are correct to 
within error estimates. The result for the subtracted integral is roughly proportional to M 2 , 
but for M 2 near 0.2/j, 2 or larger, terms even beyond M 4 appear important. 

The range of A 2 values used in the fits was from 155/i 2 to 200yU 2 in steps of 5/x 2 . For 
A 2 < 150/x 2 there is some distortion. For A 2 < 120yU 2 there is significant distortion, largely 
due to the ^ integral, which is badly approximated by the few points that satisfy the cutoff. 

The number of Fock states required for Pauli-Villars particles is approximately 1.5 times 
the number for physical states. A listing of counts for two cases is given in Table |T[ Making 
/ii larger does decrease the number of Pauli-Villars states but this increases the coefficients «j 
and thereby amplifies errors in the integrals. Also, with fewer states, the integrals themselves 
are approximated less accurately. 
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20 



15 



10 






0.000 0.002 0.004 0.006 0.008 0.010 

(^/A) 2 

FIG. 2. Subtracted one-loop fermion self energy. The Pauli-Villars masses are n\ = 10/i 2 , 
/x| = 50/i 2 , and = 100/i 2 . The solid lines are from an analytic expansion in M 2 given in ( |2. 10 ) 
of the text; additional terms are needed for M 2 > 0.2/i 2 . 
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TABLE II. Number of Fock states used in two typical cases. 



Pauli-Villars 



A 2 


K 


N ± 


physical 


= 10/u 2 


t4 = 50^ 2 


/i 2 = 100^ 2 


total 


200 n 2 


20 


25 


25975 


22602 


11142 


3305 


37049 


200 Li 2 


24 


30 


44943 


39162 


19293 


5695 


64150 



III. A SOLUBLE MODEL 

We now turn to the consideration of a nonperturbative problem. 



A. The effective Hamiltonian 



An effective Hamiltonian of the sort investigated by Greenberg and Schweber [[| and by 
Glazek and Perry |16] can be obtained from the Yukawa Hamiltonian [14|] by modifying the 
momentum dependence in the fermion kinetic energy, (M 2 + p 2 L )/p + — > (Mq + M^p + ) / P + , 
and by keeping only the no-flip three-point vertex in a modified form where the longitudinal 
momentum dependence is simplified. The fermion kinetic term in the Hamiltonian has a 



structure similar to that of the self-induced inertia term shown in Eq. (C.2) of Ref. |14 
This is a generalization of a static source. We include one Pauli-Villars field, which will 
prove sufficient in this case. The resulting light-cone Hamiltonian = P + P~ E is given by 

dp + d 2 p± 



H 



cff 
LC 



167T 3 p 4 



-{Ml + M'p + )Y,blb p _ a 



+ P 
+ 9 



dq + d 2 q± 
16n 3 q + 
dpfd 2 p^ 1 



1Qtt 3 Pi 



H 2 + <ll t , Mi + gj t 

; a q a q I T a lq a lq 



q-r 2. _ g-r 

dp\d 2 px_2 f dq + d 2 q± 



167T 3 P2" 



167r 3 g" 







pi) 


a q 
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(pi^ 




\ptj 





-2 2- 



q) 



with p = (p + ,p±) and 



16vr V% - ([) , 
{h<r^ py ,} =16ir 3 p + 8(p-p?)5, 
The Fock-state expansion of an eigenvector is 

dqfd 2 q_[_ 



n,ni 



^16n s p 



in 

n 



ni 



'^(p-p-J2ii-J2- 



- i ' \ l(i~ : 'V/,' j I 

{n ' ni \q v r_^p) — 



dr~jd 2 r^j 



167r 3 rj 



y/n\n\\ 



n ni 

&iRn«4i°> 



(3.1) 



(3.2) 



(3.3) 
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The normalization condition for this state is 

• $ a = 1Qtt 3 P + 5(P - P) , 
which yields the following condition on the individual amplitudes: 



n „ m „ 

1 = E II / d^dVi II / dr}d 2 r ± 



n,ni i 



(n,ni) 



(^•^-E^-E^ 



(3.4) 



(3.5) 



We seek a solution to 



B. Analytic solution 



With i/i = qf / P + and Zj = rj~/P + , the amplitudes must then satisfy 

{n ' ni \q v r v p) 



^E 



1 / 



P / 

+ \ T 



167r 3 g+ VP + + 4 
dr + d 2 r± ( p + — r + " 



-M\/ni + 1 



(n,ni+l) 



(g i ,L j ,r,p-r) 



+ 



E 



+ \ T 



(3.6) 



(3.7) 



167r 3 r+ \P + ?i 



{n,ni 1 (q r 1: . . . ,rj_i,Tj +1 , . . . ,r ni ,p + rj) > . 



By construction, this coupled set of integral equations is identical in basic form to the 
equations considered by Greenberg and Schweber ||. Their factorized ansatz for a solution 
suggests that we try 



dn,ni) 



z (-9) n (-ig) ni (p + Ytt v_i 



n 



(3.8) 



This is indeed a solution, provided that M = M and 

r/ _ <7 2 /P+ln Ml //x 
16vr 2 7 + 1/2 ' 



(3.9) 



Although 7 can be assigned any of a range of values, 1/2 is the natural choice, and we will 
use this value for the remainder of the paper. With this choice, the one-boson amplitude is 



proportional to Jy(l — y). 
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The normalization condition ( [3.5|) implies 



E 



Z fi£ (2n + 2m + l)!n!m! (16tt 2 )"+™i 



(3.10) 



Thus we can fix the bare mass and the wave function renormalization. However, there 



remains the bare coupling.^ 



C. Coupling renormalization 

To fix the coupling we use (:0 2 (O):) = $£. :0 2 (O): <3> CT . For the analytic solution this 
expectation value reduces to 



(:0 2 (O):> = E 



2Zn 



(g/^) 2n (g/^ 



\ (2n + 2n 1 )\n\n 1 \ (16tt 2 )"+"i 



(3.11) 



From a numerical solution it can be computed fairly efficiently in a sum similar to the 
normalization sum 



- ni . / n 

(:0 2 (O):) = En=i,n 1 =onr / ^Axifl / dr/dV^- £ 



(3.12) 



x 



(n,ni) 



(g^n^Z-E^-E^ 



With the bare parameters determined, we "predict" a value for the slope of the fermion 
no-flip form factor. It is related to the transverse size of the dressed fermion. From [I7| we 
find a useful expression for the form factor 



F(Q 2 ) 



1 



2P^ 



(P + M |J + (o)|Pt) 



£eW 167^(1- 5>W£fcu)II 



(3.13) 

lg7r 3 ^i^.Pl^nf^.Pl.). 



dxid 2 pj_ 



where the matrix element has been evaluated in the frame with 

M 2 

P=(P+,p- = — ,0 ± ), Py =(0,p- = 2p^P/P + , Py± ), Q 2 =p 2 ± , 
ej is the charge of the jth constituent, and 



(3.14) 



Pn 



P±i - XiP~f± i ^ J 
P±i + (1 - Xi)p~,± i = j ■ 



(3.15) 



2 In this model the bare coupling is finite. 
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A sum over Fock states is understood. 

When the fermion is assigned a charge of 1, and the bosons remain neutral, the analytic 
solution yields 



F(Q 2 ) = ZZ n , ni 



(g 2 /167T 3 ) n+ni /-i 



„ i n ni 

/*(i-I>-£ 



n\n\\ 



(3.16) 



n 



yidyid 2 q ±i 



n 



Zjdzjd 2 r±j 



with 



q'±_ = q±~ yp~<± and r' ± = r±- zp^± . 
The slope is extracted as 

Z(n/^ 2 + nJtt) (g/^g/^ 



no) = -£ 



^ (2n + 2m + 3)!n!m! (IBtt 2 )^ 1 
Numerically, one can compute F'(0) from 

^'(0) = E II / ^V*II / drffrjj 



(3.17) 



(3.18) 



(3.19) 



n,ni j 



X 



^fvL + E^vi,U^)( 2i) r i; z-E^-E^ 



x^^r^P-E^-E- 



with represented by finite differences. It turns out that single derivatives of typical 
amplitudes can be better approximated than the double derivatives in the Laplacian. Inte- 
gration by parts in ( |3.19| ) then leads to a computationally better quantity 



j 



F> (V) = ~ E II / dqt^q±iU I dr+d 2 r ±] 
n,ni i 

E 

% 

+E 

3 

which differs from F'(0) by surface terms which vanish as A — > oo. 



(3.20) 



fv^^^r^p-Ei-E^) 
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D. Distribution functions 



To further explore the wave functions (fr( n ' ni \ we compute distribution functions for the 
constituent bosons 



f B (y) = En, m IT / dqftfqu n • 1 / dr+d 2 r ±j E % - <?i7^ H 



x 



(n,ni) 



B1) (2,£,;£-Ei«-£; 



and the Pauli-Villars boson 



(3.21) 



3=1 



X 



' (n,ni) 



(g^i^z-E^-E^- 



Their integrals yield the average multiplicities 

(n B ) = I fB(y)dy, (n PV ) = I f PV (z)dz. 



For the analytic solution ( |3.8|) we obtain 



fs(y) 



and 



(n B ) 



fMv) = E 



^ (2ra + 2n x - l)!n!ni! (16vr 2 ) n+n i 



(npv) = E 



(<7//^(<7/^i) 2ni 



fi ) ^ (2n + 2ni + l)!n!ni! (16vr 2 ) n+n i 

For a numerical solution, the integrals can be approximated by sums. 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



IV. DLCQ APPLIED TO THE SOLUBLE MODEL 
A. Discretization 

The basic momentum discretization and approximation of integrals are discussed in 



Sec. |II B 1| . From these we construct discrete approximations to the eigenvector the 
coupled equations ( |3.7p for the amplitudes, and the derived quantities (:0 2 (O):), F' and 
distribution amplitudes. Creation operators for discrete momenta are defined by 



61. 



--bt 



tt/L 



(4.1) 
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such that they satisfy simple commutation relations 



On'n 



These follow from ( |3.2|) and the discrete delta-function representation 



S(p- p ') = — f^L\ 8 n , n . 
The discrete approximation of the eigenvectors is then 



(4.2) 



(4.3) 



n ni 



L i 



= Vs^k E E II E II E W: a+ E 

n, ni n j=l j=l 



X 



i ?i ni 

(n ' ni) (m,, a) -rriL n <&. II 4 Jo) 



where 



(n,ni) 



27T / 7T x 2 ' 



(n+m)/2 



( (n,ni) 



(4.4) 



are rescaled amplitudes, for which the normalization condition ( |3.5|) becomes 

n ni 

i = E n E n E k £ - £ m, - £ k 

n,ni i=l rrj, i j=l / j j 



(4.5) 



The most convenient basis for a numerical calculation is the number basis (or oscillator 
basis), which eliminates summation over states that differ by only rearrangement of bosons 
of the same type. We define collections of sums with a prime ri7=i EL as being restricted 
to one ordering of the momenta and introduce factorials iV{™ } = N m \N m2 \ ■ ■ ■ where iV m 
is the number of times that m l appears in the collection {mj. The amplitudes for this 
number basis are 



(n,m) 



n\ri\\ 



^(n,m) 



(4.6) 



with normalization 



n m 



i = EnE'nE'r Bi) 

n,ni i=l m ; j=l / 



In this basis the discretization of the coupled equations ( |3T7|) yields 



(4.7) 
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TABLE III. Basis sizes for DLCQ calculations in the soluble model with parameters M 2 = /x 2 , 
/i 2 = 10/U 2 , and A 2 = 50/u 2 . The numbers of physical states are in parentheses. 



K 



N± 


3 


5 


7 


9 


11 


13 


15 


17 


1 


3 


8 


18 


38 


36 


65 


110 


185 




(2) 


(4) 


(7) 


( 12) 


( 19) 


(30) 


(45) 


( 67) 


2 


19 


70 


218 


265 


590 


1120 


822 


1410 




( io) 


( 32) 


( 127) 


( H9) 


( 343) 


( 754) 


( 453) 


( 626) 


3 


43 


222 


958 


1408 


4460 


17031 


22486 


21635 




( 22) 


( 102) 


( 367) 


( 736) 


( 2671) 


( 9230) 


(13213) 


(13531) 


4 


75 


872 


3714 


9259 


49394 


50966 


110254 


328966 




(38) 


( 330) 


( 1399) 


( 5913) 


(32363) 


(32124) 


(55319) 


(172247) 


5 


99 


2028 


13702 


54100 


95176 


386140 


1553576 






(50) 


( 722) 


( 5699) 


(28065) 


(66371) 


(232400) 


(1038070) 




6 


139 


3982 


35666 


126748 


536758 


2907158 








(70) 


( 1548) 


(12991) 


(69245) 


(391511) 


(2107688) 






7 


195 


7734 


79794 


519325 


1317392 










(98) 


( 2780) 


(32891) 


(276299) 


(1008539) 








8 


275 


11736 


172118 


1165832 












(138) 


( 4268) 


(61947) 


(687394) 











1 



E 



1 f!L_ JHY ^ n+1 ' ni \m h m,l_ j ,n-m) 



(4.8) 



E 



N. 



rrk\ + 



n \t 



rrii 



N {L,L} fn-V 



n 



^(n-l,ni) 

^+^( m ,l_.l,n-l) 



Mi, ■ ■ ■ ,2B,-_i, ffi*4-i, • • • ,mndj,n + rni) 



Nn.v 



{l ' 3Y (^) 7 ^ _1) (mi.ii.---^-i.W.---.U.a + if)| 



where n = K_ — J2i Hk ~ J2jLj, {zHi}' is the set of boson momenta without and a tilde 
implies division by /i except for Lj_ = fiLj_/ir. This is a matrix eigenvalue problem, which 
for given g, fi, M, /ii, and A we solve for ip and Mq. The cutoff A 2 is applied as a limit on 
the invariant mass of individual particles, rather than on the total invariant mass of a Fock 
state. Typical basis sizes are given in Table pi . 

The bare coupling g is fixed by setting a value for 



(:0 2 (O):> -EIIE'IIE'E^ 

n,m i=l nil j=l I k=l k 



(4.9) 



The combination of the matrix equation and the imposed constraint on (:</> 2 (0):) is solved 
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iteratively. The form factor slope F'(0), the distribution functions, and average multiplicities 
are all approximated by similar discrete sums over the amplitudes ?p( n ' ni \ 



B. Numerical techniques 



The matrix equation (|4.8|) is solved using the Lanczos algorithm |18| for complex sym- 
metric matrices, which is a special case of the biorthogonal Lanczos algorithm |T9| , pO| . Given 
an initial guess U\ for an eigenvector of a complex symmetric matrix A, a sequence of vectors 
{u n } is generated by the following steps: 

v n+1 = Au n - b n u n _ x (with 61 = 0) 

O-n — v n+l " u n 

V 'n+1 = V n+1 ~ a n U n (4.10) 
bn+1 = \fv' n+1 ■ V' n+1 

u n +i = v' n+1 /b n+1 . 

The dot products do not involve conjugation, and the constants a n and b n are in general 
complex. The process will fail if b n+ \ is zero for nonzero v' n+l , which can happen in principle 
but does not seem to happen in practice [pOR . If v' n+1 is zero, the process terminates naturally. 



The vectors u n _i, v n+ \, v' n+1} and u n+ i can all be stored in the same array. At any one 
time only two vectors, one of these and u n , need to be kept. 

The vectors u n are orthogonal to each other, and the a n and b n form the diagonal and 
co-diagonal of a complex symmetric tridiagonal matrix which represents A in the basis {u n }. 
If the process has terminated with a v' n+1 = for some n, the tridiagonal representation is an 
exact representation for some subspace, and diagonalization yields some of the eigenvalues 
of A. If the process is terminated at some arbitrary early point, the eigenvalues of the 
tridiagonal matrix will approximate those of A. The approximation is particularly good 
for the extreme eigenvalues after only a few iterations. Depending on the initial guess, the 
number of iterations may need to be only 20, independent of the size of A. To reconstruct 
the eigenvectors of the original matrix, all of the u n need to be kept. Because only two are 
needed in the Lanczos algorithm, the others can be written temporarily to disk and retrieved 
later. 

We use the analytic solution ( |3.8| ) as the initial guess. Its components are either real 
or imaginary, and the process of matrix multiplication and division/multiplication by b n 
preserves this structure in a controlled way. The diagonal elements be shown to be 

real and the off-diagonal elements b n are either real or imaginary. This reduces the storage 
needed and eliminates the need for explicit complex arithmetic. 

To further reduce storage requirements, we take full advantage of the transposition sym- 
metry and sparsity of the matrix. Only nonzero elements and their indices are stored. The 
coupling g is factored out so that the matrix can be reused without change in the iterations 
that solve for g. 

To improve convergence, we include weighting factors of the sort discussed in Sec. |11 B 1| 
and Appendix ||. The circular form for transverse sums is used for two-body amplitudes and 
the extended trapezoidal rule is used for all others, with one exception. If the coefficients 
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( |B2|) for the extended trapezoidal rule become negative, a rectangle approximation is used. 
This is because restoration of symmetry for the weighted matrix requires the square roots 
of the weights. Schematically the symmetrization process is 

X] '/ "'./"./ = > y/WiWjAij^/uTjUj = £y/WiUi , (4-11) 

j j 

where ^WjWjAij is the new symmetric matrix. 

The complete specification of the weighting factors requires selection of integration order, 
because the limits of integration are interrelated by the cutoff. The simplest reduction 
of these interrelationships is made if all sums in one transverse direction are done before 
those in the orthogonal transverse direction, and all of these before the longitudinal sums. 
Within each of these three groupings the sums are done for one particle at a time in the 
ordered momentum list. One consequence of this choice is that the transverse directions are 
treated asymmetrically (except for the two-body sectors). This induces a small transverse 
asymmetry in the amplitudes of the solution, including the two-body amplitudes. The 
asymmetry disappears in the numerical limit N±_ — > oo. 



C. Results 

We have solved the discrete eigenvalue problem (|4.&Q for various cases. The physical 



parameter values chosen were 7 = 1/2, M 2 = fi 2 , fx\ = 10/i 2 , and (:0 2 (O):) = 1 or 2. 
The parameters that control the numerical approximation, namely A 2 , K, and N± (or L±), 
were varied to study convergence with basis sizes up to ~ 520, 000. The ranges of these 
numerical parameters are shown in Tables |IV| and [V|. The transverse scale L± was chosen 



such that N± radial points satisfy the invariant-mass cutoff for one-boson states at the value 
of longitudinal momentum that yields maximum transverse width. The bare fermion mass 
M was allowed to vary from its analytic, infinite-cutoff value of M in order that M could 
be held fixed. The tables list the values of M along with those of the bare coupling g, as 
set by Q4.9p, the average boson multiplicity (n B ), and, for (:</> 2 (0):) = 1, the slope F'(0) of 
the fermion form factor. The analytic, infinite-cutoff values are also included. 

The values of the form factor slope are very poor approximations. This is due to the 
sensitivity to N± of the finite difference representation of the derivatives in (|3.20j ). A good 



approximation requires at least A^i = 8, which implies very large basis sizes even for small 
K. 

The results for g and (%) are surprisingly insensitive to variation in K and N_\_. Only 
the cutoff A 2 is important. This is mirrored in the distribution functions shown in Figs. 
Again, variations in K and N± make little difference; however, one can see that the cutoff 
has an important effect. Smaller cutoffs produce an enhancement in the interval (0.4, 0.8).[| 

The amplitude for the one-boson state is shown in Fig. [7[ The analytic result is shown 
for comparison. As can be seen, the two shapes are nearly identical. 



Recall that the distribution function does not have a fixed normalization but instead determines 
the average multiplicity, which is then also enhanced at finite cutoff. 
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TABLE IV. Numerical parameter values and results from solving the 


model ei£ 


renvalue prob- 


lem. 


The physical parameter values were M 2 = 


/j? for the fermion mass, n\ = 


10/r 2 for the 




-Villars mass, 




5> 2 (0):) = 1 to fix the coupling g. 






(a/m) 


2 K 


_l_ 


/iL ± /n 


(M //i) 


g/n 


(fin) 
Vol 


100/z 2 -F"(0) 


50 


11 


4 


0.8165 


0.8547 


13.293 


0.177 


-0.751 


50 


13 


4 


0.8165 


0.8518 


13.230 


0.172 


-l.Uld 


50 


15 


4 


U.olOO 


U.o4lUo 


i q etc 


0.178 


-0.715 


50 


17 


4 


0.8165 


0.8289 


13.392 


0.180 


-0.565 


50 


y 


5 


1.2062 


0.8601 


14.023 


0.179 




OU 


9 


6 


1.2247 


0.8377 


14.323 


0.179 


-0.582 


50 


9 


7 


1.4289 


0.8302 


14.386 


0.179 


-0.658 


50 


9 


5 


1.2062 


0.8601 


14.023 


0.179 


-0.547 


100 


9 


5 


0.7143 


1.0520 


12.565 


0.174 


-0.239 


200 


9 


5 


0.5025 


1.1980 


10.191 


0.172 


-0.139 


oo 




analytic 


1.0000 


13.148 


0.160 


-0.786 






TABLE V. Same 


as Table 


IV 


except (:c/> 2 (0):) = 2. 


(A/M) 


2 K 




N ± 


LlL ± /tt 




(Mo//.) 2 






50 


11 




4 


0.8165 




0.5068 


21.541 


0.368 


50 


13 




4 


0.8165 




0.5166 


21.327 


0.352 


50 


15 




4 


0.8165 




0.4496 


22.323 


0.366 


50 


17 




4 


0.8165 




0.4439 


21.930 


0.364 


50 


9 




5 


1.2062 




0.5340 


22.396 


0.367 


50 


9 




6 


1.2247 




0.5109 


22.507 


0.369 


50 


9 




7 


1.4289 




0.5204 


22.287 


0.366 


50 


9 




5 


1.2062 




0.5340 


22.396 


0.367 


100 


9 




5 


0.7143 




0.9353 


20.962 


0.359 


200 


9 




5 


0.5025 




1.3080 


18.034 


0.347 


oo 






analytic 






1.0000 


19.420 


0.308 
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FIG. 3. The boson distribution function fs at various numerical resolutions, with (:</> 2 (0):) = 1 
and A 2 = 50/U 2 . The solid line is the analytic result. 
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FIG. 4. The boson distribution function fs for different cutoff values, with (:0 2 (O):) = 1 and 
numerical resolution set at K = 9 and N± = 5. The solid line is the analytic result. 
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FIG. 5. Same as Fig. | but for (:0 2 (O):) = 2. 
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FIG. 7. The one-boson amplitude t/X 1 ' ) as a function of longitudinal momentum fraction y 
and one transverse momentum component q x in the q y = plane. The analytic result is shown 
in (a) and the numerical result in (b) for A 2 = 50 fj 2 , K = 17, and N± = 4. Both correspond to 
<:0 2 (O):> = 1. 



24 



V. CONCLUSION 



In this paper we present a new method for the renormalization of Hamiltonian light- 
cone-quantized field theories that maintains Lorentz invariance and other symmetries. The 
main difficulty which is confronted by such methods is the construction of the counterterms. 
We employ the traditional generalized Pauli-Villars method ||. With a sufficient number of 
Pauli-Villars fields, perturbation theory is regulated while Lorentz symmetries and discrete 
symmetries are preserved with a minimal number of counterterms. These counterterms are 
generated automatically. We hypothesize that these counterterms are sufficient to regulate 
the nonperturbative problem. 

In Yukawa theory, Pauli-Villars regularization preserves chiral symmetry [p,L^| , unlike 



the invariant-mass regulator. A similar outcome arises in QED where one needs Pauli-Villars 
regularization of the vacuum polarization loop to recover its vanishing at q 2 — > plj. These 



examples show that a covariant method is necessary in the nonperturbative context to find 
all counterterms. 

Given that the theory is finite by suitable Pauli-Villars regularization, we can impose a 
regulator to limit the Fock space so as to produce a tractable numerical problem. A finite 
matrix approximation can then be obtained with use of the DLCQ procedure |l|]. In the 
finite matrix problem we face the numerical difficulties of non-Hermitian matrices and large 
basis sizes. These difficulties are successfully addressed in a 3 + 1-dimensional model || 
constructed to have an analytic solution. This model requires one Pauli-Villars boson as a 
regulator. We also study the DLCQ approximation to the one-loop fermion self energy in 
Yukawa theory, where three Pauli-Villars bosons are needed . 



The non-Hermitian matrices are handled by the complex symmetric Lanczos diagonal- 



ization algorithm fiq-fZOf . This technique is ideal for the extraction of extreme eigenvalues 
and their eigenvectors. It takes full advantage of the sparsity of the Hamiltonian matrix. 
For a given basis size, storage requirements are minimized. 

The basis sizes required in the calculation are reasonable. The presence of Pauli-Villars 
particles, at the chosen mass and cutoff values, increased the model basis by only 100% 
and the loop-calculation basis by 150%. Given the sparsity of the matrix, increases of 
these magnitudes are quite acceptable. However, smooth convergence and extrapolation 
from bases of minimal size require the introduction of special integral weighting methods to 
DLCQ. The dramatic improvement which can occur is illustrated in Fig. p]. 

With these methods we have obtained agreement between the numerical and analytic so- 
lutions of our model. The convergence of the numerical result in longitudinal and transverse 
resolution is remarkably rapid. The result is sensitive only to the cutoff used to limit the 
Fock space, but even there the convergence to the analytic result is clear. The methods seem 
well suited to situations where low-mass states have a small mean number of constituents. 



The natural next step is to extend the model toward a more realistic theory, namely 
Yukawa theory. The fermion can be given proper dynamics, and Yukawa-type interactions 
can be reintroduced. Once Yukawa theory itself can be studied with our nonperturbative 
method, there may be useful applications to the Higgs sector of the Standard Model. We 
are sufficiently encouraged by the success of the Pauli-Villars program for the examples 
discussed here to believe that it will have general applicability to QCD in 3 + 1 dimensions. 
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APPENDIX A: LIGHT-CONE COORDINATES 



We define light-cone coordinates by 

x ± = x° ± x 3 , 



(Al) 



with the transverse coordinates x± = (x 1 , x 2 ) unchanged. Covariant four-vectors are written 
as e.g. x^ = (x + ,x~, x±), with the spacetime metric 



Explicitly, 



[0 2 \ 

2 

0-10 

\0 -I J 



x-y = g^y u = -{x + y +x y + )-x ± - y ± . 



We also make use of an underscore notation: for position-space variables we write 

x = (x~,x ± ), 

while for momentum-space variables 

k = (k + ,k ± ). 

Then 



k ■ x = -k + x — k± ■ x I . 
2 



Spatial derivatives are defined by 

a 4 



9_ 







0; 







dx + dx~ dx l 

The gamma matrices = 7 ± 7 s = (7 T )^ satisfy the familiar relation 

{7 M ,7"} = 2^ 

with g^ u the light-cone metric . It is simple to verify that the (Hermitian) matrices 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 
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satisfy 



A| = Ad 



A+A 



0. 



A,+A_ = 1, 



(A9) 



(A10) 



so that they serve as projectors on spinor space. In the Dirac representation of the 7- 
matrices: 



A, 



which has two eigenvectors, both with eigenvalue +1: 





(I 





1 


\ 


1 





1 





-1 


2 


1 





1 







\0 


-1 





1 / 



(All) 



1 

V2 





1 

voy 



X- 



1 

V2 



( \ 
1 



v-iy 



(A12) 



These serve as a convenient spinor basis for the expansion of the field ip + = A + ip on the 
light cone. 



APPENDIX B: WEIGHTING METHODS 

New weighting schemes have now been developed for use with the DLCQ grid. They are 
based on extensions of the trapezoidal rule and Simpson's rule to the situation where the 
integration domain does not end on a grid point; they are also related to open Newton-Cotes 
formulas. The basic approach is to derive formulas for one-dimensional integrals and then 



iterate them for higher-dimensional integrals [15 



The extended trapezoidal rule is obtained from consideration of an integral from, say, 
Xq to £3. The relevant graph is shown in Fig. The grid points are at X\ and x 2 , which 
are separated by a standard spacing h. The other points are at the integration domain 
boundaries at distances of hi and from the grid points. The integral of a function / is 
then approximated by 



I-X3 

/ fdx~a 1 f(x 1 )+a 2 f(x 2 ), (Bl) 

Jx 



with 



ax = (h + h L + h R )(h + h L - h R )/2h , 

a2 = (h + h L + h R )lh + h R -h L )/2h. (B2) 

The coefficients dj are chosen to provide exact results for linear functions. The standard 
trapezoidal rule is recovered when h L = h R = 0. If h L = h R = h, a standard open 
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FIG. 8. Spacing of grid points for an arbitrary function. 

Newton-Cotes formula results. When the extended rule is combined with the standard rule 
for interior intervals, a general composite rule is obtained. The extended rule is then used 
twice, once at each end, with Hr or Hl set to zero. 

The extended Simpson's rule follows from similar steps. Two forms are needed, one for 
three grid points and another for four. Any situation with more grid points can be handled 
with a composite rule obtained by combining these rules with the standard Simpson's rule. 
For the three-point case, consider an integral from x to £4, with grid points at x±, x 2 , and 
x 3 . The regular grid spacing is h; the extra points at the beginning and end are separated 
by hi and Hr, respectively. The approximation to the integral is then 

/ fdx ~ a l f{xi) + a 2 f(x 2 ) + a 3 f(x 3 ) , (B3) 

where 

fll = (4/i 3 + 12h 2 h L + 9hh 2 L + 2h\ + 3hh 2 R + 2h R ) /12h 2 , 

a 2 = (Ah 3 - 3hh 2 L -h\- 3hh 2 R - h%) /3h 2 , (B4) 

a 3 = (4h 3 + \2h 2 h R + 9hh 2 R + 2h R + 3hh 2 L + 2h\) /12h 2 . 

These coefficients are constructed to provide exact results for quadratic functions. When 
h-L = h>R = they reduce to the coefficients found in Simpson's rule, and, because of the 
greater symmetry, the rule becomes exact for cubic functions as well. 
The four-point rule is also exact for cubic functions. It is 

/ fdx ~ axf(xi) +a 2 f(x 2 ) +a 3 f(x 3 ) + a 4 f(x 4 ) , (B5) 

Jx 

where 
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a x = (9h 4 + 24h 3 h L + Tlh 2 h\ + 8hh 3 L + h 4 L - Ah 2 h 2 R - 4hh 3 R - h R ) /2Ah 3 , 

a 2 = (27h 4 - 36h 2 h 2 L - 20hh 3 L - 3h\ + lSh 2 h 2 R + 16hh 3 R + 3h R ) /24h 3 , (B6) 

a 3 = (27h 4 - 36h 2 h 2 R - 20hh 3 R - 3h R + l%h 2 h 2 L + 16hh 3 L + 3h A L )/2Ah 3 , 

a 4 = (9h 4 + 2Ah 3 h R + 22h 2 h 2 R + 8hh 3 R + h R - 4h 2 h 2 L - Ahh\ - h\) /2Ah 3 . 

These integration formulas greatly reduce the size of the errors, as shown for the extended 
trapezoidal rule ( |B1|) in Fig. but they do not result in systematic behavior. The lack of 
systematic dependence is primarily due to the use of a square grid to approximate a circular 
domain in the transverse direction. One way of putting this is that the iterated Cartesian 
integrals try to approximate ir (a key factor in the area of the circle) as well as approximate 
the integral itself. 

To overcome this square/circle problem, the integral is written in polar coordinates 

J dxdyf{x, y)= l - J* d(r 2 )f(r 2 , 0) . (B7) 

The points of the square grid lie on circles of varying radii shown in Fig. |1(]. The r« are 
easily computed from the coordinates of the square grid. The squares of these radii are used 
as the grid points for a trapezoidal approximation to the radial (r 2 ) integral. Because the 
limit R 2 does not fall on one of these points, the extended trapezoidal rule (|BT| ) must be 
used for the last interval. Clearly, the intervals are not of equal length; however, they are 
on average of order 3h 2 , where h is the spacing in the square grid. For the first 10 circles, 
the average spacing in r 2 is actually closer to 2h 2 . 

The number of points on the square grid that fall on any one circle come in multiples of 
4, because of reflection symmetries. These points can be used to approximate the angular 
integral on each circle via another application of the trapezoidal rule. For the self-energy 
integral J, however, the integrand is independent of angle and one can simply use one point 
or average the values at all points. The contribution to the weighting of a grid point is then 
the same for all transverse points on the same circle. 

The circular weighting in the transverse direction can be combined with either the ex- 
tended trapezoidal rule or the extended Simpson's rule in the longitudinal direction. A 
comparison of the two is shown in Fig. |TT[ The relatively large excursions for small iV^ are 
due to the small number of grid points involved for this case of a large boson mass. For 
larger iV^ the extended Simpson's rule is seen to result in less excursion between different 
values of K, and is preferred for the self-energy integral. Results for the extended Simpson's 
rule are compared with those of the ordinary DLCQ sum in Fig. [I]; in this case the results 
for the extended trapezoidal rule would not be visibly worse. 
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FIG. 9. Same as Fig. |l] except that the numerical values oscillating about the analytic answer 
are computed without transverse circular weighting and with only trapezoidal weighting in the 
longitudinal and transverse directions. The lack of circular weighting destroys the smoothness of 
the results shown in Fig. |]. 
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0.00 0.05 0.10 0.15 0.20 



1/N ± 

FIG. 11. One-loop fermion self energy. Results for the extended trapezoidal and Simpson 
rules are compared. The horizontal line is the exact result. The DLCQ grid parameters take the 
ranges K = 20, 21, . . . , 25 and N± = 5, 6, . . . , 30. Points calculated with the same value of K are 
connected by lines. 
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